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ABSTRACT 

Aims. An interesting question of contemporary cosmology concerns the relation between the spatial distribution of galaxies 
and dark matter, which is thought to be the driving force behind the structure formation in the Universe. In this paper, we 
measure this relation, parameterised by the linear stochastic bias parameters, for a range of spatial scales using the data of the 
Garching-Bonn Deep Survey (GaBoDS). 
Methods. The weak gravitational lensing effect is used to infer matter density fluctuations within the field-of-view of the survey 
fields. This information is employed for a statistical comparison of the galaxy distribution to the total matter distribution. The 
result of this comparison is expressed by means of the linear bias factor b, the ratio of density fluctuations, and the correlation 
factor r between density fluctuations. The total galaxy sample is divided into three sub-samples using i?-band magnitudes and 
the weak lensing analysis is applied separately for each sub-sample. Together with the photometric redshifts from the related 
COMBO-17 survey we estimate the typical mean redshifts of these samples with z — 0.35, 0.47, 0.61, respectively. 
Results. Using a fiat ACDM model with Om — 0.3, Qa ~ 0.7 as fiducial cosmology, we obtain values for the galaxy bias on 
scales between 1' < 6ap < 20'. At 10', the median redshifts of the samples correspond roughly to a typical comoving scale of 
3, 5, 7 h~^Mpc with h = 0.7, respectively. We find evidence for a scale-dependence of b. Averaging the measurements of the bias 
over the range 2' < 6a,p < 19' yields b = 0.81 ± 0.11, 0.79 ± 0.11, 0.81 ± 0.11 (Icr), respectively. Galaxies are thus less clustered 
^ ] than the total matter on that particular range of scales (anti-biased). As for the correlation factor r we see no scale-dependence 
• • , within the statistical uncertainties; the average over the same range is f = 0.61 ± 0.16, 0.64 ± 0.18, 0.58 ± 0.19 (Icr), respectively. 
^ ' This implies a possible decorrelation between galaxy and dark matter distribution. An evolution of galaxy bias with redshift is 
not found, the upper limits are: Ab < 0.2 and Ar < 0.4(1(t). 
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1. Introduction 

In comparison to the total mass in the Universe, galax- 
ies take - considering their mass - only a minor part in 
the big picture of structure formation. They formed from 
the baryonic component, embedded in the fluctuations of 
the dark matter density field, whose total mean density is 
very much lower than that of dark matter. Due to their 
relatively easy observability, it would be very convenient if 

* Based on observations made with ESO Telescopes at the 
La Silla Observatory 

** Founded by merging of the Institut fur Astrophysik 
und Extraterrestrische Forschung, the Sternwarte, and the 
Radioastronomisches Institut der Universitat Bonn. 



galaxies were perfect tracers - unbiased tracers - of the to- 
tal mass distribution; all statistical properties of the mass 
structure could then be derived from galaxy catalogues. 

Indeed, it is rather unlikely that galaxies are unbiased 
tracers, because the laws determining the galaxy distri- 
bution are very complex and highly non-linear. The pri- 
mordial gas from which they form requires special con- 
ditions to be able to cool and fragment into galaxies 
(White & Frenk 1991; White & Rees 1978). Due to shock 
heating of the baryons and energy feedback between galax- 
ies and the intergalactic medium, the properties of the 
gas feeding galaxy formation gradually changed with time. 
Furthermore, galaxies interact with each other or with the 
baryonic intergalactic medium, merge or get accreted into 
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more massive galaxies (Lacey & Cole 1993). These mech- 
anisms probably produced the large diversity in galaxy 
masses, colours, morphologies and chemistry wc observe 
today. Based on our current knowledge it would be very 
surprising if this complexity would eventually result in a 
simple, linear, one-to-one relationship between the galaxy 
density and total matter density, making galaxies unbiased 
tracers. But by studying this dark matter-galaxy relation- 
ship we can learn more about galaxies. 

Observing the relation between the invisible dark mat- 
ter field and the galaxies is a particularly tough problem. 
However, with gravitational lensing at hand, we now have 
a technique to directly unravel this relationship. The im- 
portance of "cosmic shear" as a tool for cosmology was 
proposed in the early 1990s by Blandford et al. (1991), 
Miralda-Escude (1991) and Kaiser (1992). Since these pi- 
oneering days of gravitational lensing the techniques and 
surveys have been refined to make valuable contributions 
to the ongoing research on structure formation on cosmo- 
logical scales. In particular, the investigation of the rela- 
tion between galaxy and dark matter distribution using 
the weak gravitational lensing effect has become almost 
standard (Scljak et al. 2005; Klcinheinrich et al. 2005; 
Sheldon et al. 2004; Hoekstra et al. 2003; Guzik & Seljak 
2001; McKay et al. 2001; Fisher et al. 2000; Hudson et al. 
1998; Braincrd et al. 1996). This paper will focus on the 
lensing technique as well. 

1.1. Quart tifying ga laxy bias 

From the point of view of statistics, quantifying galaxy 
bias leads to the question how one can parametrise dif- 
ferences in the statistical properties - not the obvious 
differences between two particular realisations ~ of two 
random fields. Both the distribution of galaxies and the 
distribution of dark matter are thought to be realisations 
of statistically homogeneous and isotropic random fields. 
Commonly, one uses a parametric way to describe the bi- 
asing between two (random) density fields, for instance 
galaxy and matter distribution or the distributions of two 
different galaxy populations. Biasing between two density 
fields, say pg and pm, can in general be quantified using the 
joint probability distribution function (PDF) P ((5g, i5ni) of 
the density contrasts (density fluctuations) 



at the same point in the density fields at some redshift (lo- 
cal Eulerian bias). The probability of finding density con- 
trasts of 6g and within an interval of d^g and dS^, re- 
spectively, equals P{Sg,S^)d5^ddg. The density contrasts 
are smoothed to a certain scale, R, before investigating 
their PDF. Often the special type of smoothing kernel is 
set by the method that is used to determine the bias. To 
lowest, second order (first order moments vanish by defi- 
nition, (5g) = (^m) = 0) the only relevant parameters for 



biasing are the scale-dependent parameters of the "linear 
stochastic bias" 

m - M ; riR) = -^M±= , (2) 

V ^^^^ ^{SiHSD 

which differ from unity in the case of two biased fields. 
These two parameters provide a complete description for 
the bias between Sg and 5yn for Gaussian fields only and 
therefore are a full description only on large smoothing 
scales where linear perturbation theory applies (on "linear 
scales" ) . The linear bias factor 6 is a measure for the dif- 
ference in clustering strength. The correlation coefficient 
r, on the other hand, measures partly the stochasticity 
in the relation between the density contrasts, for instance 
how well minima and maxima of the density fields co- 
incide. Only partly, because the correlation coefficient is 
also sensitive to the non-linearity in this relation, which 
becomes relevant on smaller smoothing scales where the 
fields are usually no longer Gaussian. In fact, they cannot 
be Gaussian since ^ > — 1 by definition. Despite this de- 
generacy, the linear stochastic bias parameters are clearly 
defined on non-linear scales; it is just their interpretation 
which is no longer straightforward. To disentangle a non- 
linear but deterministic relation between two density con- 
trasts from a stochastic one requires higher-order statis- 
tics, like for example in the context of the "non-linear 
stochastic bias" scheme (Dekel & Lahav 1999). 

1.2. Galaxy bias in observations 

Observationally, galaxy bias can be derived from the one- 
dimensional PDF of galaxies, P{Sg), (Marinoni et al. 
2005; Sigad et al. 2000), redshift space-distortions (Pen 
1998; Kaiser 1987), weak gravitational lensing (Seljak et 
al. 2005; Sheldon et al. 2004; Pen et al. 2003; Hoekstra 
et al. 2002; Wilson et al. 2001; van Waerbeke 1998; 
Schneider 1998) and counts-in-cells statistics (Conway et 
al. 2005; Tegmark & Bromley 1999; Efstathiou et al. 1990). 
Additionally, the large-scale fiow of galaxies can be used to 
make a POTENT reconstruction of the total mass field on 
large scales which can be compared to the galaxy distribu- 
tion (Sigad et al. 1998; Dekel et al. 1993). Gravitational 
lensing (Schneider et al. 2006; van Waerbeke & Mellier 
2003; Bartelmann & Schneider 2001) provides a promising 
new method in this respect becaiise it allows for the first 
time to map the total matter content (mainly dark mat- 
ter) independent from the galaxy distribution. The work 
of this paper is based on this technique. 

A brief overview of the current status of the obser- 
vational results is given in the following. Note that the 
given conclusions to some extent depend on the assumed 
cosmological model. We quote only the conclusions for the 
concordance ACDM model (cf. Tegmark et al. 2004). In 
the local universe, galaxies are almost unbiased tracers 
on linear scales of about 8/i.~^Mpc and larger (Seljak et al. 
2005; Verde et al. 2002; Lahav et al. 2002; Loveday et al. 
1996). However, this is probably not true on smaller scales. 
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A comparison of the theoretical dark matter clustering - 
which is constrained by the cosmic microwave background 
anisotropics, gravitational lensing and the Lyman-a forest 
(Tegmark et al. 2004) - and the observable galaxy clus- 
tering suggests that on smaller scales ~ l/i~^Mpc galaxies 
are less clustered than the dark matter ( "anti-biased" ) be- 
coming positively biased, & > 1, on even smaller scales be- 
low ~ O.l/i-^Mpc. Hoekstra et al. (2001, 2002) use in their 
work on the VIRMOS-DESCART survey (van Waerbeke 
et al. 2001) and RCS (Gladders & Yee 2001) weak grav- 
itational lensing to measure the linear stochastic bias for 
galaxies with a median redshift of — 0.35, covering a 
range from 0.1/i^^Mpc to 6.3ft.^^Mpc. They claim to have 
observed such a dip in the linear bias factor. Also based 
on gravitational lensing there is evidence that the ratio 
hjr stays close to unity from sub-megaparsec scales up 
to - S/i-^Mpc (Sheldon et al. 2004; Hoekstra et al. 2002; 
Guzik & Seljak 2001; Fisher et al. 2000), thus from non- 
linear to linear scales. The analysis of the bispectrum of 
the galaxy clustering in the 2dFGRS (CoUess et al. 2001) 
led Verde et al. (2002) to the conclusion that on scales 
between 5/i~^Mpc and 30ft.~^Mpc the biasing relation be- 
tween dark matter and galaxies is essentially linear (see 
also Lahav et al. 2002). The same conclusion was drawn 
several years earlier by Gaztanaga & Frieman (1994) based 
on the APM survey (Maddox et al. 1990). However, re- 
cently the work of Kayo et al. (2004) has questioned a 
strict linear relation on scales < 10/i^^Mpc by studying 
the three-point correlation of galaxy clustering as a func- 
tion of morphology, colour and luminosity, this time in the 
SDSS (York et al. 2000). 

Subdividing the galaxies into various subsets gives a 
more detailed picture of galaxy biasing. At low redshift, 
clustering is a function of morphological type, spectral 
type, colour and luminosity (e.g. Madgwick et al. 2003; 
Zehavi et al. 2002; Norberg et al. 2001; Benoist et al. 1996; 
Tucker et al. 1997; Loveday et al. 1995; Davis & Geller 
1976). Late- type, blue, spiral or star forming galaxies are 
less clustered than early-type, elliptical or red galaxies 
with a relative linear bias factor of about &iod/&biuc ~ 1-4 
on scales of roughly 8 /i~^Mpc (e.g. Wild et al. 2005; 
Conway et al. 2005; Zehavi et al. 2002; Norberg et al. 
2002; Baker et al. 1998). On large scales, the relative bi- 
asing between red and blue galaxies does not seem to be 
well described by a simple linear biasing function which 
according to Wild et al. (2005) (see also Conway et al. 
2005) is ruled out with high significance using counts-in- 
cells statistics in redshift space. Wild et al. observe a scale- 
dependent non-linear bias between red and blue galaxies 
with a dominant stochasticity component for typical phys- 
ical scales of about 7/i~^Mpc up to 31ft.~^Mpc. Relative 
bias between red and blue galaxies is therefore both non- 
linear and stochastic. It is therefore also expected that at 
least for some galaxy populations the bias with respect 
to the dark matter distribution is non-linear and stochas- 
tic as suggested by simulations (Yoshikawa et al. 2001). 
This, however, has not been measured directly so far. 
Observational evidence for the relation between red and 



blue galaxies being non-deterministic was already given 
some years ago by the work of Tegmark & Bromley (1999) , 
which was based on the Las Campas Redshift Survey 
(Shectman et al. 1996), and by Blanton (2000). Moreover, 
galaxy bias seems to be a function of redshift (Marinoni 
et al. 2005; Magliocchetti et al. 2000) which is expected 
both from simulations (Weinberg et al. 2004 and refer- 
ences therein) and analytical models (Tegmark & Peebles 
1998; Mo & White 1996; Fry 1996). 

In this paper, we apply the method of Hoekstra et al. 
(2002; Sect.EI) to the Garching-Bonn Deep Survey (Sect. 

to obtain the linear stochastic bias coeSicients h and r 
of three galaxy subsets binned by their apparent i?-band 
magnitude. The selection of the galaxy samples is outlined 
in Sect. 121 After presenting our results in Sect. [S] we close 
with a discussion and conclusions in Sect.El We will start 
with a brief introduction to the formalism of weak gravi- 
tational lensing and the aperture statistics here employed. 

Unless otherwise stated we use a ACDM model with 
a„ = 0.3, = 0.7 and iJo = /i 100 kms"! Mpc"^ with 
h — 0.7. A scale-invariant in = 1, Harrison-Zel'dovich) 
spectrum of primordial fluctuations is assumed. As trans- 
fer function, encoding the physical properties of the cold 
dark matter fluid, we use Bardeen et al. (1986). 

2. Formalism 

2.1. Weak gravitational lensing 

Weak gravitational lensing uses the shapes of distant 
galaxies - the source galaxies or, as we will also call them, 
background galaxies ~ to infer the distribution of the to- 
tal matter. This is based on the fact that light is deflected 
by density fluctuations so that the tidal gravitational field 
of the matter density inhomogeneities along the line-of- 
side towards a galaxy changes the shape of its image. We 
consider only cases in which the hght rays emitted from 
a source galaxy traverse only regions in space with rela- 
tively small perturbations in the gravitational fleld (weak 
lensing regime); this holds for almost every galaxy. 

Theory. In the weak lensing regime the differential distor- 
tion effect of the tidal gravitational fleld is well described 
by a two-dimensional linear mapping over the whole ap- 
parent size of one galaxy. The relevant components of the 
linear transformation A are the convergence k, which mag- 
nifies or demagnifies the size of a galaxy, and the shear, 
7 = 7i -|- i72 , which stretches the image of a galaxy along 
some direction: 

m = (^-^-^' -\ ) , (3) 

V ^^"^ 1 - K + 71 y ' ^ ' 

where B is an angular position on the sky. Note that for 
convenience 2D vectors are written as complex numbers, 
as for instance B = Q\ ^iQ^- 

According to the theory of weak gravitational lensing, 
convergence and shear are, to lowest order, a projection 
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of the three-dimensional density contrast (5m of the matter 
in the Universe via 



dwWiw)SUfKiw)e,w) 



(4) 
(5) 



where 



Aw pb [w ) — „ , — (6) 



2c^a (w) /k [w') 

is a weight given to the density contrast at a comoving ra- 
dial distance w from the observer. The functions a(w) and 
/k{w) are the cosmic scale factor and the comoving an- 
gular distance, respectively, at comoving distance w. The 
variable is the comoving horizon size. By Phiw) we 
denote the distribution of the background sources in co- 
moving distance. The shear 7 as a function of the direction 
on the sky 6 is called the cosmic shear map. According to 
Eq. (0) it is a convolution of the convergence k, thus the 
projected matter density contrast, with some kernel. 

In the formalism presented here, we assume that we 
always observe only small patches of the celestial sphere; 
small enough to approximate the topology of the patch 
by a tangential, Cartesian plane (flat-sky approxima- 
tion). This is a very good approximation for the sur- 
vey fields considered here, which are smaller than Ideg^. 
Furthermore, the 3D coordinate system, such as for S^i in 
Eq. 10}, is chosen such that w is a comoving distance along 
some fixed reference line-of-sight and f}^(w)9 a 2D-vector 
perpendicular to the reference line-of-sight; w is also used 
as look-back time parameter to account for the fact that 
6m is a function of time. 

Connection to the real world. To quantify the shape of 
a galaxy one defines the complex ellipticity, e = ei -I- ie2, 
which is related to the quadrupole moments of the 
light distribution in the galaxy image (e.g. Schneider 
et al. 2006). It transforms under rotations according to 
e' = exp {—2iip) e where if) is the rotation angle. 

Seitz & Schneider (1997), for example, showed that 
under the action of the linear transformation A the in- 
trinsic ellipticity, i.e. the unlensed galaxy ellipticity, Egj of 
a source galaxy is transformed into the image ellipticity, 
e, according to 



Cs +7 



(7) 



if 7 ^ 1 which is fulfilled in the weak lensing regime. The 
practical importance of this equation stems from the fun- 
damental assumption that galaxies have intrinsically no 
preferred direction and are therefore randomly oriented, 
i.e. (eg) = 0. This makes the observed ellipticities of galax- 
ies, e(0), unbiased estimators of the cosmic shear field in 
the direction 0: 



{m) = i{e) 



(8) 



In this picture, the shear is a function of the conver- 
gence, and the convergence is related to the projected (5m. 



Consequently, it should be possible to make a reconstruc- 
tion of the total matter distribution based on the observed 
ellipticities of source galaxies, or a reconstruction of struc- 
tural parameters such as spatial correlation functions. 

Of course, galaxies are in general not intrinsically 
round objects, eg 7^ 0; the ellipticities of single galaxies 
have typically (|es|) ~ 0.3 with a comparable scatter of 
cr(|es|) ~ 0.3. This makes them in fact very noisy estima- 
tors of the shear, considering that the shear signal induced 
by gravitational lensing is typically about one percent of 
this value. Therefore, the average over many galaxy ellip- 
ticities is required in weak lensing applications. 

2.2. Aperture statistics 

Aperture mass. For our application here, which aims at 
the linear stochastic bias parameters on some smoothing 
scale, we are not interested in the convergence field itself 
but in the convergence field smoothed to some typical scale 
using an aperture filter u. A convenient quantity for this 
purpose is the aperture mass Afap (Schneider et al. 1998) 
defined as 



Map(0ap,e) 



1 

ap 



W -0\ 



(e') 



(9) 



The variable 0ap is the aperture radius setting the filter 
scale. It has been shown that if the filter u is compensated, 
J dx X u (x) = 0, then the aperture mass can be deter- 
mined from the shear field itself (Schneider 1996, 1998) 



A^ap/x (0. 



e-e'\ 



7t/x (0') , (10) 
2 /"^ 

q{x) = — \ dssu{s)-u{x) , (11) 
Jo 



where 

7t {e') 

7x {0') 



-Re (7(6>')e-2'"^ 
-Ini {j{e')e-^'"^ 



(12) 
(13) 



denote the tangential and cross shear component with re- 
spect to the aperture centre 6, respectively. The variable 
4> is used as definition for the argument of the complex 
number ■3 = 9' -0, i.e. '&i + iz?2 = li^le'"^- 

In this definition, the actual aperture mass or so-called 
E-mode of the aperture mass is obtained by using the 
tangential shear (with respect to the aperture centre), 7t, 
while choosing the cross shear, 7x , gives the B-mode, , 
of the aperture mass. As the shear originates from a sin- 
gle scalar field, k, the two shear components are related 
to each other (cf. Schneider et al. 2002). Therefore, not 
all conceivable shear field configurations are produced by 
gravitational lensing. The allowed configurations of 7 are 
called E-modes, while the other independent configura- 
tions are called B-modes. For that reason, a signature of 
B-modes is used in this paper as an indicator for system- 
atics in the data reduction, especially the point-spread 
function (PSF) correction, which has to be performed to 
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compensate the instrumental and atmospheric influence 
on the galaxy image. Note, however, that on scales smaller 
than about a few arcmin a non-zero B-mode can be pro- 
duced by intrinsic alignments of the source galaxies (e.g. 
Heymans et al. 2004; Hirata et al. 2004) or spatial clus- 
tering of the source galaxies (Schneider et al. 2002). 

Since the ellipticity of a galaxy at is, in the weak 
lensing regime, an unbiased estimator of the shear 7, one 
could construct an estimator for the aperture mass that 
can be directly applied to a galaxy catalogue in order to 
obtain a Map-map for some survey field (Hoekstra et al. 
2001). We are, however, interested in the relation between 
matter and galaxy distribution in a statistical sense. As 
we will see soon, for this purpose it is not even necessary 
to make an actual map ~ even though this could be a 
possible strategy. Before we discuss aperture statistics we 
introduce a quantity to analyse the spatial distribution of 
galaxies. 

Aperture number count. In a similar fashion to the aper- 
ture mass, we can define (Schneider 1998; van Waerbeke 
1998) the aperture number count N{6, ^ap) which mea- 
sures the fluctuations of the galaxy number density with 
the same filter u as A/ap for the convergence field: 

A^(0,0ap) ^ ^ j A^O' u (^^^) 5n{e') (14) 

where n{6) and n denote the (projected) number density 
of galaxies in some direction and the mean number den- 
sity of galaxies, respectively. The quantity 5n = n/n — 1 
is the projected number density contrast of the galaxies. 

Usually, the galaxies probed with N{0, 6'ap) and those 
galaxies used to construct Mi^p{6,9ap) maps are differ- 
ent; the latter tend to be more distant as they probe the 
matter field containing the galaxies used for N. For that 
reason, we call the "N-galaxies" foreground galaxies and 
the "Map-galaxies" background galaxies. 

The projected galaxy number density contrast is re- 
lated to the three-dimensional galaxy number density con- 
trast, Sg, via 

5n{e)^ dwpf{w)6g{fK{w)e,w) . (16) 
Jo 

The function pt (w) is the distribution of foreground galax- 
ies in comoving distance w which will be estimated from 
the observed distribution in (photometric) redshift. Note 
that the galaxy distribution in redshift, pf{z), and distri- 
bution in comoving distance are related by: 

Pfiw) = pKz)^ = Pfi^w)) ' , (17) 

where z{w) is the redshift as function of w. 

Eq. (|16|1 is the counterpart to Eq. It is the pro- 
jected density contrast of the galaxy density, while k is 
the projected density contrast of the total matter. 



2.3. Aperture statistics and correlation functions 

In order to estimate the linear stochastic bias, Eq. 
using the aperture number count {N related to dg) and 
aperture mass statistics (Map related to 6^) we need to 
estimate the second-order moments of the aperture statis- 
tics, i.e. (Ar"(6'ap)M™(6'ap)) with m + n = 2. There are 
two principal ways to estimate the 2'^'^-order moments of 
N and Map: either by placing apertures at different po- 
sitions onto the field (Hoekstra et al. 2001), or indirectly 
by estimating and transforming the two-point correlation 
function of the galaxy number density, cosmic shear and 
their cross-correlation (Hoekstra et al. 2002). In this pa- 
per, we are going to use the latter method. How the cor- 
relation functions relate to the aperture statistics will be 
summarised in the following. 

Power spectra. From the statistical point of view, the 
joint 2"'^-ordcr moments of the aperture statistics are 
fluctuations, {N^ i&ap)) and (M^ (^ap)), and correlations, 
{N{9ap)Map{0ap)), of smoothcd (statistically homoge- 
neous and isotropic) random fields. They are therefore 
auto- and cross-correlation power spectra seen through a 
filter [I{x)]'^ (Hoekstra et al. 2002): 

poo 

{Mlp ((?ap)) = 27r / d«P. {t) [I {m,p)f , (18) 

{N (0ap) Map (0ap)) = ^TT / d« Pn« {£) [I (^^ap)]' , (19) 

Jo 

poo 

{N^d^p)) = 27r / dUPn{£)[I{ee,p)]\ (20) 
Jo 

with the filter function 



I{x)= / dssu{s)Jo{sx) , (21) 
where 

(27r)2fe(£ + £')P.(|^|) = {m~<i')), (22) 

(2^)2fc(£ + /)P,„(|^|) = {5fi{£)k{£')) , (23) 

{2TT)H^{£ + £')Pnm = {SmSh{£')) , (24) 



are the convergence auto-correlation, P^, convergence- 
galaxy number density contrast cross-correlation, P^n, and 
the galaxy number contrast auto-correlation power spec- 
trum, Pn. Su denotes the Dirac delta function and a tilde 
is used indicates the Fourier transform, such as 

k{£) = j d^eK{e)c+''^^ . (25) 

We use Jn{x) for the n*^-order Bessel function of the first 
kind. 

Using Limber's equation in Fourier space (Kaiser 1992) 
we can derive these power spectra from Eq. H16|l and Eq. 
O: 

P,{t) ^ rdu,&ll!p„/— ^, J , (26) 
Jo [/k(w)]2 V/k(w) J 
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[/k(w)]2 



/k(w) 



/k(«^) ' 



(27) 
(28) 



They will be needed below for the calculation of the cali- 
bration factors. Here, Pin(fc,w) represents the 3D matter 
power spectrum as a function of comoving distance w. In 
the above equations, the (3D) galaxy-matter cross-power 
spectrum, 

(27r)3jD(fc + fc')-Pm.„(fc,u.) = {K{k,w)~5^{k' ,w)) , (29) 
and (3D) galaxy power spectrum, 

(27r)3jD(fc + fc')-Pn(fc, w) = (4(fc, w)4(fc', , (30) 

are expressed in terms of the linear stochastic bias in 
Fourier space: 



1,2/7 \ Pn\k,w) . , 



fm,n(fc, W) 



VPn(fc,H-Pm(fc,U') 



(31) 



They are the Fourier space counterparts of the linear 
stochastic bias parameters in Eq. The Fourier space 
bias parameters are uniquely defined - they are uniquely 
attached to Fourier modes which represent a certain phys- 
ical scale -, whereas the real-space definition Q requires 
further specification of a (smoothing) window function, 
Wr(x), through which the density fluctuations are ob- 
served. Once the Fourier space bias is known and a window 
function is defined we can always work out the real-space 
bias parameters: 



and 



/dfcfc^ (b^F.n)(fc,C.^W^R(fc)P 
/dfcPP,„(fc,C^)|W^R(fc)|2 

1 /dfcP(6rP.n)(fc,oyR,(fc)p 
h{R,uj) /dA:fc2p^(A:,w)|Vl?R(fc)|2 



(32) 



(33) 



where WR(fc) is the Fourier transform of the window 
Wyi{x) with respect to x. Hence, all information on the 
linear stochastic bias (2"'^-order biasing) is contained in 
the Fourier space bias parameters. 

Unbiased galaxies have r = 6 = 1 for all and w. 
In general, however, they are time- and scale-dependent. 
Note that for a very narrow, or even delta function like 
Pf(u') Eq. H28|l diverges because the assumptions made 
in the Limber equation break down for those cases (see 
Simon 2006). 

Correlation functions. The 2"'^-order statistics of random 
fields is completely described by the power spectra. 
Equivalently, we can consider correlations between pairs of 
points in the random fields in real space which gives rise to 
the two-point correlation functions. Although not carrying 
different information than the power spectra, correlation 



functions are of great practical value because they are rel- 
atively easily calculated; we just have to consider pairs of 
galaxies in our case. 

These correlators are here: a) the angular correlation of 
the foreground galaxy positions, uj{9), b) the mean tan- 
gential shear about foreground galaxies, (7t) (0), and c) 
the shear-shear correlations ^± (6*) as determined from the 
ellipticities of the background galaxies (Hoekstra et al. 
2002): 



(7t) (e) 



{Sn {9 + x) 5n {x)) 

°° ds s 

^P„ {s) Jo {s9) , 
27r 

{Sn {9 + x) 7t (x)) 

°° ds s 

Pk„ {s)J2{s9), 



2tt 



{-/t{9 + x)-tt{x))±{^^ {9 



(34) 
(35) 

+ x)-i^{x)) (36) 



The correlator lij{9) is a measure for the probability of 
finding a galaxy at a separation 9 from another galaxy. 
The galaxy-galaxy lensing, (7t)(^), is the mean tangential 
shear around a foreground galaxy at a separation 0, and 
the two-point shear-shear correlations, S,±{9), quantify the 
correlations of the cross- and tangential shear components 
relative to the line connecting two background galaxies 
with an angular separation 9. All correlators are linearly 
related to the corresponding power spectra. 

Transformation integrals. The relations between the corre- 
lators and the power spectra, Eqs. (|M|l -lp!H )l . can be in- 
verted with respect to the power spectra (Schneider et al. 
2002; Hoekstra et al. 2002), so that invoking Eqs. ((TH|l - 
(|20|l we can express the aperture moments in terms of the 
two-point correlators: 



(^4 (^ap)) 



da; x [^+ (a;6'ap) P+ [x) (37) 



+ ^_(a;0ap)r-(a;)], 



{N (0ap) A4p (0ap)) 



{N^ (^ap)) 



Ax X (7t) (a;6'ap) F (x) , 



Ax X uj {x9a_p) r+ (x) , 



where we introduced the auxiliary functions 



r+,. (x) EE (27r)^ 



As s [I (s)] Jo, 4 (sx) , 



(38) 
(39) 

(40) 
(41) 



The Eqs. H37|) - (|39|l are the basis for the method that is 
used in this paper. The two-point correlators are estimated 
from the data and then afterwards integrated in order to 
obtain the 2"'*-order moments of the joint PDF of aperture 
mass and aperture number count. 



F{x) = (2'Kf / Ass [I{s)fj2{sx) 
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Fig. 1. Frequency of apparent galaxy i?-band magnitudes 
in the shallow, medium and deep part of GaBoDS (fore- 
ground galaxy samples). The distribution functions have 
been normalised by the area between 18 < i? < 24 mag 
since galaxies fainter than 24 mag are not considered in 
this paper. As can be seen in this plot, all three data sets 
have roughly comparable distributions for R < 24 mag. 

If we substitute in the foregoing Eq. (|37|) and (|38|l tan- 
gential shear components, 7t, by cross shear components, 
7x , and vice versa then we obtain the transformation in- 
tegrals for the corresponding B-modes; we then have (7x) 
in (|38|l and the ^_ term in H37|) changing sign. 

3. Data 

Here we will give only a brief account of the GaBoDS. For 
details concerning the GaBoDS, its data reduction and 
catalogue creation, we refer the reader to Schirmer (2004), 
Erben et al. (2005) and in particular Hetterscheidt et al. 
(2006). 

3.1. The GaBoDS fields and their reduction 

The GaBoDS comprises roughly 18.6 deg^ of high-quality 
data (seeing better than one arcsec) in i?-band taken 
with the Wide Field Imager (WFI) mounted on the 2.2m 
telescope of MPG/ESO at La Silla, Chile; the 33' x 34' 
field-of-view is covered with 8 CCD chips. Due to the 
dither pattern applied, the effective field-of-view can be 
as large as roughly 40' x 40'. The data set was compiled 
mostly from archival ESO data, for which the archive util- 
ity querator (Pierfederici 2001) has been developed, to- 
gether with about four square degree coming from our 
own observations. The positions of the fields were chosen 
randomly from regions of small stellar densities at high 
galactic latitudes. The limiting magnitudes of the fields is 
inhomogeneous, ranging between 25.0 mag and 26.5 mag 
(5(7 in a 2" aperture radius) in the iZ-band depending on 
the exposure time and on the fraction of time the see- 
ing was acceptable for gravitational lensing applications. 
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Fig. 2. Magnitude vs. half-light radius plot of objects 
found by SExtractor in one particular field. Stars ap- 
pear as almost vertical branch and can be separated from 
galaxies with high confidence. The solid and dashed box 
roughly encircles objects excluded for the lensing cata- 
logue (Schirmer et al. 2003). 

The data set can roughly be categorised into a shallow 
{t < 7ksec, total 9.6 deg^), medium (7ksec < < < lOksec, 
total 7.4 deg^) and deep (lOksec < i < 56 ksec, total 
2.6 deg^) set depending on the total usable integration 
time t for each field. 

The data imposed new, high demands on the data re- 
duction, which resulted into the development of a data 
reduction pipeline whose usage is not restricted to the 
aforementioned instrument only; it has successfully been 
tested on data from various other instruments (Erben et 
al. 2005). 

For the final analysis, we consider the shallow, medium 
and deep part of the GaBoDS comprising in total 52 WFI 
fields corresponding to an area of about 15 deg^. We re- 
jected nine fields: CAPO and all fields belonging to the 
CO series (eight), as the quality of the PSF correction was 
decided to be not sufficient enough for weak gravitational 
lensing applications (see Hetterscheidt et al. 2006). As will 
be explained shortly we do not consider galaxies, for both 
lensing and foreground object catalogues, that are fainter 
in the i?-band than 24 mag. Applying this cut at the faint 
end makes the GaBoDS categories shallow, medium and 
deep roughly comparable with each other as can be seen 
by the magnitude histogram in Fig. ^ and it allows us to 
estimate the redshift distribution of galaxies (see below). 

3.2. Selection of lensing catalogues 

After the data reduction process, SExtractor 
(Bertin & Arnouts 1996) was used to compile a catalogue 
of source galaxy candidates needed for the cosmic shear 
analysis. For the rather conservative selection of source 
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candidates, the final co-added science frames are first 
smoothed with a Gaussian kernel of 2.5 pixel FHWM. 
One pixel corresponds to 0'.'238. A source candidate fur- 
ther needs to consist of at least 5 contiguous pixels with a 
total flux greater than 1.5 a above the background noise 
level, and it has to possess a clearly defined quadrupole 
moment (cl ^ 0, analyseldac) and centroid. Stars and 
galaxies are distinguished in a magnitude vs. half-light 
radius plot of the selected objects (see Fig. |21 for an 
example). In this scatter plot, stars that are not too faint 
are clearly identified as a column of objects with roughly 
identical half-light radius r^,. Objects with a half-light 
radius smaller than are rejected as source candidates. 
An exception are objects in the faint part (fainter than 
23.5 mag in i?-band) near this column. 

As accurate measurements of galaxy shapes are the 
key in a weak lensing analysis, the quadrupole moments 
in the galaxy light profiles of the source candidates have to 
be corrected for PSF effects; atmospheric turbulence and 
instrumental effects also distort the galaxy images. This 
is done using the KSB method (Kaiser et al. 1995). A 
detailed description of the PSF correction procedure may 
be found in Erben et al. (2001), Heymans et al. (2006) 
or Hetterscheidt et al. (2006). The PSF fitting polynomial 
used is of order two or three. 

In the estimators of the aperture statistics, every 
source galaxy is weighted with a statistical weight. This 
weight, Wi, is defined by the variance in cUipticity 
of the 12 nearest neighbours of a galaxy i in the mag- 
nitude vs. half-hght radius diagram: wi = 1 / {a"^ + a^) , 
where — 0.16 is the variance of the unlensed galaxies. 
In the case that the PSF corrected ellipticity of a galaxy 
exceeds |e| = 1.0 it automatically is attributed the weight 
zero and is hence not considered further in the analysis. 
Applying this cut removes rare outliers with unrealistic el- 
lipticities, produced by the KSB technique. The final lens- 
ing catalogue is split into three magnitude bins BACK, 
BACK-II and BACK-III, see Table [H 

3.3. Selection of foreground objects 

The actual foreground objects, of which the bias parame- 
ters are measured, are selected with the same SExtractor 
parameters as the galaxies candidates in the lensing cata- 
logue. Galaxies are finally selected from this catalogue via 
a manually defined box in the magnitude/half- light radius 
diagram around the star branch. 

In order to select for the bias analysis different mean 
redshifts of the (foreground) object catalogues, we subdi- 
vide the object catalogue into the three different i?-band 
bins FORE-I, FORE-II and FORE-III as stated in Table 

m 

3.4. Distribution in redshift of the galaxy samples 

To estimate the redshift distribution of the galaxies - both 
foreground objects and background sources - we average 



Table 1. The table lists the limits of the magnitude bins, 
the total number of objects for all 52 fields (deep, medium 
and shallow fields in GaBoDS), the mean redshift and the 
Icr-variance inside each bin. 

Foreground object catalogue 

sample bin limits [mag] ^objects (2) 

FORE-I 19.5 <R< 21.0 6.5 x 10* 0.34 ± 0.18 
FORE-II 21.0 <R< 22.0 1.2 x 10^ 0.47 ± 0.22 
FORE-III 22.0 <R< 23.0 2.5 x 10^ 0.62 ± 0.27 

Background source catalogue 

sample bin limits [mag] ^objects (z) 

BACK 21.5 <i?< 24.0 6.2 x 10^ 0.67 ± 0.29 

BACK-II 22.0 <R< 24.0 5.5 x 10^ 0.69 ± 0.28 
BACK-III 23.0 <R< 24.0 3.5 x 10^ 0.74 ± 0.28 



the photometric redshift distribution in the different mag- 
nitude bins of the fields A901, AXAF and Sll (see Fig. 
O . These three fields are contained in the deep part of the 
GaBoDS and were observed as part of the COMBO-17 
survey (Wolf et al. 2004) in 17 colours yielding quite 
accurate photometric redshifts with an uncertainty of 
Sz « 0.02 (1 -I- z) for objects brighter than i? = 23 mag. 
Less accurate but still available are photometric redshifts 
for objects with 23 < i? < 24. The photometric redshift 
distribution of the COMBO-17 galaxies is assumed to be 
representative for our whole catalogue. 

High-redshift galaxies with z > 1.4 are "missing" in 
the COMBO-17 sample because they were reassigned a 
redshift z < 1.4. Recently, in Coe et al. (2006), the Hubble 
Ultra-Deep Field has been used to validate the photomet- 
ric redshifts in COMBO-17. It has been found that the 
agreement is good for R < 23.7 and especially tight for 
R < 23. We conclude, therefore, that we have got a reli- 
able estimate of the redshift distribution in our samples. 

For the source galaxies carrying the A/ap-signal, the 
magnitude bin BACK is used throughout. As can be seen 
in Tabled by varying only the lower limit, but keeping 
the upper limit of the magnitude bin fixed to i? = 24 mag, 
one cannot shift the mean of the background redshift dis- 
tribution to much higher values than z « 0.7; essentially, 
only the number of sources in the bin decreases. A large 
mean redshift of the source galaxies is desired to achieve 
a good lensing efficiency but more important is a large 
number of galaxies to achieve a good signal-to-noise ratio. 
Since we do not use objects fainter than 24 mag in order 
to maintain good accuracy in the estimate for the red- 
shift distribution of the background and to have a roughly 
homogeneous data set, the bin BACK for all three fore- 
ground bins FORE-I, FORE-II and FORE-III is the best 
choice. 
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Fig. 3. Redshift distribution of the foreground and background galaxies as estimated from the photometric redshifts 
in the COMBO-17 fields A901, AXAF (CDFS) and Sll (dashed doted lines); the histograms arc not normahsed to 
unity. The solid lines are maximum- likelihood fits of Eq. H42(l to the histograms. 



The COMBO-17 sample used to estimate the redshift 
distribution in the galaxy sub-samples is relatively small. 
Clearly, it has features - large galaxy clusters or voids - 
which are not representative for the whole GaBoDS sam- 
ple. For example, consider the peaks at low redshift in the 
foreground samples. Fig. |31 In order to have a smoother, 
more representative distribution we fit an empirical red- 
shift distribution to the COMBO-17 histograms: 



1 



— cxp 



(42) 



The best-fits are shown in Fig. |3| They are used in the 
analysis further on instead of the COMBO-17 histograms. 
The best fit parameters are compiled in Table El These 
parameters are highly degenerate for which reason their 
statistical errors are not given in the table. However, sta- 
tistical uncertainties of the mean redshift in the different 
samples and uncertainties of the galaxy bias parameters 
originating from uncertainties in p{z) will be given in the 
following. 

Clearly, the estimated redshift distributions still suf- 
fer from cosmic variance errors because the COMBO-17 
survey area is with ~ 0.75 deg^ relatively small. In order 
to get an estimate for the statistical uncertainties due to 
cosmic variance in the samples' redshift distribution we 



use the widely applied Jackknife method: The photomet- 
ric redshift distributions of merely two of the three fields 
are combined. With three ways of combining this yields 
overall = 3 Jackknife samples. To estimate the standard 
deviation of the mean redshift, z, one computes from each 
Jackknife sample the mean redshift, Zi. According to the 
Jackknife method the statistical Icr-error of the mean is 
then roughly: 



2,, A^-1,-^, ,2 



(43) 



where z is the mean redshift obtained by combining all 
three COMBO-17 redshift distributions. The results for 
a{z) are listed in Table HJ As can be seen there the un- 
certainty of z ranges from (j{z)/z « 10% to cr{z)/z « 2% 
for FORE-I to BACK, respectively. This behaviour makes 
sense because the number of galaxies increases when going 
from the shallower to the deeper samples. 

The problem of the calibration of redshift distributions 
for cosmic shear studies has recently been studied by van 
Waerbeke et al. (2006). They find a statistical uncertainty 
of (t{z) = 0.03 — 0.04 for a 0.75 deg^ survey with mean 
z ^ 1. This value is somewhat higher than our estimate. 

The Jackknife samples can also be used to assess how 
the statistical uncertainty of the full p(z)'s translates into 
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Table 2. Best- fit parameters of the template redshift dis- 
tribution, Eq. (gH), to the COMBO-17 histo grams, z is the 
mean of the template redshift distribution. The statistical 
errors of z are derived from the field-to-field variance in 
COMBO-17. 



galaxy sample 


zo 


Q 







z 




FORE-I 


0.534 


0.509 


3.173 


0.35 


± 


0.03 


FORE-II 


0.765 


0.617 


5.839 


0.47 


± 


0.03 


FORE-III 


0.945 


0.830 


5.103 


0.61 


± 


0.02 


BACK 


1.069 


0.809 


7.369 


0.68 


± 


0.02 


BACK-II 


1.072 


0.988 


7.655 


0.70 


± 


0.02 


BACK-III 


1.073 


1.611 


8.560 


0.77 


± 


0.02 



the inferred galaxy bias parameters. This problem will be 
addressed in Sect. 14.31 

4. Outline of the method 

The approach to obtain the bias parameters from lensing 
adopted here proceeds in several steps: 

1. estimating the binned correlators (7t)(^) and 
C±(^) in all individual survey fields, 

2. numerical integration of the correlators to obtain 
(iV'"(6lap)M^p(e'ap)) for TO + n = 2 (E-modes and B- 
modes), 

3. repetition of 1. and 2. with bootstrapped data sets to 
obtain statistical errors of the aperture statistics in the 
single fields, 

4. combining the individual field measurements and eval- 
uating the bias parameters as a function of aperture 
radius from the combined signal (includes calibration) , 

5. bootstrapping of the combined signal to estimate the 
error in the final signal and the covariances between 
the different bins. 

A detailed account of these steps is given in the following. 

4.1. Practical estimators for the correlators 

The correlators are estimated by using 



DD DR 

uj{e) = — - - 2— + 1 



(7t) {0) = 



i± (0) 



where 



RR RR 



Ef'' w.w,X, (9) 



1 ioT 9 < \9,-e, \ <9 + 69 
otherwise 



A,, (9)=^ ^ ^ ''^^ 



(44) 
(45) 

(46) 
(47) 



The Wi are statistical weights of the individual galaxies 
which are used to account for the fact that the values 



of the image ellipticities, e^, of the galaxies do not have 
all the same accuracy. Ellipticities of fainter and smaller 
galaxies are determined with a lower accuracy than for 
larger and brighter galaxies. A^f and A^b are the number of 
foreground and background galaxies; Ct.i/j ^^'^ 
the tangential and cross ellipticity components relative to 
the line connecting the galaxy pair 

The estimator of the spatial correlation lu{9) has been 
introduced by Landy & Szalay (1993). It requires to count 
the number of galaxy pairs with a separation between 9 
and 9 + 69, namely the number of pairs in the data (fore- 
ground galaxies), denoted by DD, the number of pairs in a 
random mock catalogue, RR, and the number of pairs that 
can be formed with one data galaxy and one mock data 
galaxy, DR. The random mock catalogue is computed by 
randomly placing the galaxies, taking into account the ge- 
ometry of the data field, i.e. by avoiding outmasked re- 
gions. We make 40 random galaxy catalogues and average 
the pair counts obtained for DR and RR. 

In an estimate of uj{9), Eq. (|44|l . there is always an 
uncertainty about the mean galaxy density n which is the 
larger the smaller the area of the field under consideration. 
This introduces a bias known as the integral constraint 
(Groth & Peebles 1977), that systematically reduces the 
angular correlation, lu{9) uj{9) — C , by a constant value 
C > 0. As pointed out by Hoekstra et al. (2002) (A^) is 
independent of the integral constraint when the aperture 
filter u is, as in our analysis, compensated because 



/•OO PCO 

I dxxCT^ {x) — C / dx xT-^- 
Jo Jo 



(48) 



Therefore the estimator bias C does not make any contri- 
bution to (A^) and does not need to be corrected. This 
nice feature makes the aperture statistics with a compen- 
sated filter a convenient tool to study galaxy clustering 
even outside a weak lensing context. We will briefly come 
back to this point in Sect. 

Concerning the estimator for mean tangential shear, 
{jt){9), all pairs of foreground and background galaxies 
within separation between 9 and 9 + 69 have to be con- 
sidered; et,i is the tangential ellipticity component of the 
background galaxy with respect to the line connecting the 
foreground and background galaxy. Similarly, for ^±(^) all 
pairs of background galaxies within some separation inter- 
val are considered. 

For the GaBoDS analysis, we bin the three correla- 
tors into 800 logarithmic bins spanning a range between 
0'.'05 < 61 < 48' (the diagonal of a single WFI field). In or- 
der to reduce the computation time for the correlations, 
a binary tree data structure as in Pen & Zhang (2003), 
Moore et al. (2001) or Jarvis et al. (2004) is used. 
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4.2. Aperture filter 



Pin{k, w) by means of 



To weight density fluctuations inside apertures we use a 
compensated polynomial filter (Hoekstra et al. 2002, 2001; 
Schneider et al. 1998) 



1 



H{1 - x) 



which by definition vanishes for a; > 1; H(x) denotes 
the Heaviside step function. The filter has the effect that 
only matter or galaxy number density fluctuations from a 
small range of angular scales contribute to the N- or Map- 
signal; it acts as a narrow-band filter for the angular modes 
with highest sensitivity to £c ^ 4.25/6'ap ~ 0.68 x 27r/6'ap. 
Apertures with radius ^ap therefore effectively probe a co- 
moving physical scale of /k(w) 6'ap/0.68, if w is the median 
comoving radial distance of the galaxy sample under ex- 
amination. 

All auxiliary functions (|4U|I - 141|) . which arc required for 
transforming the correlators, vanish outside the interval 
X £ [0,2] due to the finite support of u. This reduces the 
transformation integrals l|37|) - (|39|l to a finite integration 
range [0, 26'ap]. Therefore, with square 30' x 30' WFI fields 
we are able to estimate the aperture moments out to about 
21'. 



iV2 30 



4.3. Calibration of bias parameters 

To summarise, the aperture mass Map, Eq. (|10|l . is pro- 
portional to the (weighted) projected total matter density 
contrast Jm, whereas the aperture number count iV, Eq. 
lfn|l . is proportional to the number density contrast of the 
galaxy distribution Sg. Both aperture measures are defined 
on some scale by the filter function u and the aperture size 
^ap- This is exactly what we need to study the biasing of 
the galaxy distribution with respect to the matter distri- 
bution, as has been pointed out by Schneider (1998) and 
van Waerbeke (1998). Therefore, we can define biasing pa- 
rameters in analogy to Eqs. (O (Hoekstra et al. 2002) 



6(0ap) 
'^(^ap) 



/l (^ap) 



/2 (^ap) 



{N' (M) 

V (^4^p(M) ' 

{N (gap) Map (0ap)) 



(50) 



The three-dimensional number density that N is sensitive 
to covers in general not the same volume that is probed 
by Map {pf{w) versus W{w)). Naively identifying N with 
the galaxy number density contrast, (5g, and Afap with the 
matter density contrast, S^n therefore gives the wrong bias 
parameters. This 7V/Afap-mismatch in sensitivity to den- 
sity fiuctuations at the same radial distance thus requires 
to make a correction of the bias parameters which is done 
by means of the calibration factors fi/2- 

According to Hoekstra et al. (2002), the calibration 
factors have to be calculated based on some theoretical 



Tl l^apj 



'(M2p(0ap)) 



r=b=l 



(49) /: 



2 IC'ap; 



V (^M^ap)) 
^{N^ (0ap))(M2p(0ap)) 



(51) 



{N (0ap) Map (0ap)) 



=6=1 



where (iV"(0ap)-Af™(0ap)), n + m — 2, in these equations 
have to be evaluated by Eqs. and lf^ - (P5|) . 

specifically for the redshift distributions of foreground, 
Pf(z), and background galaxies, Phiz), in the data and 
for a fiducial cosmological model assuming that galax- 
ies are not biased with respect to the dark matter, i.e. 
b{k, w) — r{k, w) = 1. 

Importantly, it turns out (van Waerbeke 1998) that 
the calibration factors /i and /2 vary only slightly, mostly 
on scales below ^ap ^5 5', for realistic aperture radii ^ap 
within a fixed fiducial cosmological model. This is strictly 
true if the dark matter power spectrum can be described 
by a power law, or - since we are, for a fixed aperture 
radius, sensitive to only a very localised range in £-space 
due to the adopted aperture filter ~ if the power spectrum 
is approximately a power law over the selected range in 
Fourier space. 

For examples, see Fig. 01 (upper left and bottom left) 
where /1/2 are plotted for three fiducial cosmological mod- 
els assuming the redshift distribution of FORE-I and 
BACK. The calibration factors show very little depen- 
dence on 6'ap. Hence, a scale-dependence of the uncali- 
brated measurements immediately indicates a real scale- 
dependence in the bias parameter without fixing the fidu- 
cial cosmology! Moreover, it means that the calibration 
factors can be worked out for the linear or quasi-linear 
regime which is understood much better than the non- 
linear regime. Still, when calibrating our measurements 
we also take into account the dependence on scale. 

We calculated the calibration factors /1/2 for a range 
of spatially flat fiducial cosmologies, fl^ + f^A = 1, using 
the redshift distribution in our data set (right column in 
Fig. assuming constraints on erg oc f^m^'^^ from clus- 
ter abundances (White et al. 1993) and the shape pa- 
rameter r — fljuh for a negligible baryon density fib ~ 
(Efstathiou et al. 1992). The relation between as and Om 
is scaled such that as = 0.8 corresponds to ilm = 0.3. This 
value of CTg for the power spectrum normalisation is sug- 
gested by the GaBoDS data (Hetterscheidt et al. 2006). 
Note that the value of erg, like for example tig = 0.9 in- 
stead of (Tg = 0.8, has virtually no impact on /1/2 and, 
therefore, on the measured linear stochastic bias. 

Predicting the power spectra, Pk, Pku and P^, requires 
a model for the redshift evolution of the 3D power spec- 
trum Pn-i{k,w). We use the standard prescription of lin- 
ear structure growth and the Peacock and Dodds (1996) 
prescription for the evolution in the non-linear regime. A 
more recent and more accurate description of the non- 
linear power spectrum is given by Smith et al. (2003). 
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Fig. 4. Left: The two figures show the scale-dependence of the calibration factors /1/2, for sample FORE-I only, for 
three diflferent fiducial cosmologies; SCDM (dotted): = 1.0, = 0; ACDM (solid): = 0.3, f^A = 0.7; OCDM 
(dashed): fim = 0.3, SIa — 0. Right: These plots were obtained by averaging fi/2 over a range of aperture radii, 
1' < ^ap < 60', assuming different fiducial cosmologies. For all cosmologies, fim is the only free parameter. The others 
are: fl^ = 1 — fii„, F — Q^h, as — 0.41 and h = 0.7. The average values for fi/2 in this figure are divided by 

fi/2{^m = 0.3), the here adopted calibration. For pf(z), we have FORE-I (solid), FORE-II (dashed) and FORE-III 
(dotted); Phiz) is as in BACK. 



Although Smith et al. predict in general more clustering 
on smaller scales than Peacock & Dodds, we found in a 
comparison between both methods only little difference 
for fi/2- 

It becomes clear from Fig. 0| that particularly the 
interpretation of the bias factor, b (calibration /i), de- 
pends on ilm. For the final calibration of the GaBoDS 
measurements we assume as fiducial cosmological model 
= 0.3, riA = 0.7, erg = 0.8, F = 0.21, h = 0.7. 

The calibration procedure outlined here works because 
correlations between fluctuations seen in N and Afap stem- 
ming from different redshifts quickly vanish with increas- 
ing mutual redshift difference. This allows us to correct 
for a mismatch in the cosmic volumes seen by N and Afap 
- for the price of making assumptions about the fiducial 
cosmology, though. In fact, the mismatch is quantified by 
the factor I//2 which is the correlation factor of the (pro- 
jected) fluctuations seen in N and those seen in Afap as- 
suming that galaxies perfectly trace mass. If this factor is 
exactly I//2 = 1 for all scales, we will immediately know 
that N and Afap probe exactly the same cosmic volume 
giving equal weight to all radial distances. For our sam- 



ples and flducial cosmology, we find I//2 = 0.95, 0.95, 0.84 
(FORE-I to FORE-III) indicating that the mismatch is 
relatively small. 

4.4. Uncertainties in the calibration factors 

Concerning the statistical errors of the calibrated galaxy 
bias parameters it has to be considered that the red- 
shift distributions of the galaxy samples are not ex- 
actly known (Sect. I3.4|l . Since the calibration just in- 
volves a rescaling by /1/2, any relative error in /x/2 re- 
sults in an equal relative error in 6(0ap) and ^(^ap). To 
estimate the relative error in /1/2, f(/i/2)//i/2i we use 
again the three Jackknife samples of the redshift distri- 
butions that have already been used in Sect. 13.41 Every 
histogram of photometric redshifts of every Jackknife sam- 
ple is fitted by the template distribution H42I) . Then, 
for our fiducial cosmology and for every galaxy sam- 
ple, we compute /1/2 for each Jackknife p(z)-template. 
Thus, for any sample we obtain a triplet of values for 
the bias calibration. Then, as in Sect. 13.41 the variance 
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between the different /1/2's is used as estimate of the Icr- 
error of the caUbration parameters. Following this pro- 
cedure, we find that the statistical uncertainty in the 
calibration due to cosmic variance uncertainties in p{z) 
is (FORE-I to FORE-III): cr(/i)//i = 16%, 4%, 3% and 
c^(/2)//2 = 4%, 2%, 8%, respectively. Thus, the relative er- 
ror is roughly cr{fi/2)/ fi/2 ^ 7%, except for /i of FORE-I 
where we have cr(/i)//i = 16%. We expect similar relative 
errors for other cosmological models. 

Another issue is the impact of uncertainties in the 
fiducial cosmological model on the calibration parameters, 
which in return influences the inferred galaxy bias. Using 
Fig. ^ we estimate the relative change in fi/2 with ~ ±7% 
(/i) and - ±2% (/a) when changing = 0.3 by ±10%. 
Therefore, having wrong by about ±10% introduces a 
systematic error into the bias parameters which is about 
±7% for b and about ±2% for r. 

4.5. Redshift and scale resolution 

The bias parameters obtained by Eqs. H50|l are in general 
averages of the true bias parameters b{k,w) and r{k,w), 
Eqs. H31|l . namely averaged over some scale k and comov- 
ing distance (redshift) w. In other words, the method ap- 
plied here has a limited resolution in redshift and scale. 
This is due to two reasons: 1. the sample of foreground 
galaxies is usually not peaked at one particular redshift 
but smeared out over some range, and 2. Icnsing is, with 
varying response, sensitive to the whole matter distribu- 
tion between a source galaxy and the observer. According 
to Hoekstra et al. (2002), the observed weighted averages, 
(6^) and (r), when using the polynomial aperture filter, 
Eq. (|49|l . are approximately 



J^'" dwhi{w,9a^p) 
/p""" Aw h3{w,9i^p) 



(52) 



(53) 



where the following weight functions have been introduced 



^3(W, 6'ap) 

with 

-Pfiltor(w,6'ap) = 

Jo 



Pfjw) 
/k(w) 
P[{w)W{w) 



-Pfiltor(w,6'ap) , 



(54) 
(55) 



fK{w)9, 



w [iie9.p)Y 



(56) 



Note that in general the central redshift and the width of 
the weights depends on the aperture radius Sap- 

We calculated the functions hi{z,9ap) and ^13(2, Sap) 
for accessible aperture radii for the redshift distributions 
in our galaxy samples and for the adopted fiducial cos- 
mological model. As typical redshift at which b and r are 
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Fig. 6. The aperture number count dispersions, as mea- 
sured in GaBoDS, for FORE-I (filled boxes), FORE-II 
(open stars) and FORE-III (open crosses). The la error 
bars have the size of the data points. Upper panel: com- 
parison to ACDM predictions assuming unbiased galaxies^ 
upper to lower fine: FORE-I (sohd), FORE-II (dashed) 
and FORE-III (dotted). Lower panel: power laws give 
excellent descriptions of the measurements. The dotted- 
dashed line denotes {N'^) as measured by Hoekstra et al. 
(2002) {A^ = 0.115, S = 0.7). 



measured we take - for every aperture radius Sap - the 
mean, Zc, of the kernels hi and /13, and as a measure for 
the redshift range over which the bias parameters are av- 
eraged we take the width, CTz, of the kernels hi and /13 
(see Fig.lSJ. With Zc being the typical redshift, the typical 
spatial scale corresponds to kc = g — f'^{w{z )) Fourier 
mode) or Rc = 27i-/fcc ~ 1.48 Sap/K(w(2c)); the redshift 
resolution achieved is roughly which also adds an un- 
certainty to the effective scale probed, aji, with a relative 
error of about cru/Rc = (Tz/zq- 



4.6. Combining measurements from different fields 

Bootstrapping I. The whole data set taken into account for 
the analysis consists of Np = 52 fields. For each field, we 
compute the two-point correlation estimates as in Sect. 
14.11 and transform them according to the integrals Eqs. 
(I37|l - (|39|l to the 2'^'^-order aperture statistics. For every 
individual field, j, we make 200 bootstrap samples of the 
foreground and background object catalogues to estimate 
the statistical error of the measurements 



a:p)e{(iV2(Sap,.)), (A^(eap.)Map(Sap,0), (Mfp (Sap..)) } , (57) 
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Fig. 5. The bias parameters estimated in this paper are averages over some redshift range, plotted here as a function of 
aperture radius; plotted is the mean, Zc, (solid Hue) and variance, Cz, (shaded area) of the weight functions hi/^{z, 6'ap) 
for each foreground sample. Mean and variance do not change significantly with aperture radius, 0ap- The average of 
Zc and Uz over all radii, Zc ± fz, is indicated in each panel. Left column: weight function for the bias factor, hi. 
Right column: weight function for the correlation parameter, ^13. 



for the various aperture radii 6'ap,i . The aperture moments 
x^^^ from all fields are combined to one final result, Xi, by 
the weighted average 



(58) 



As weight we use the reciprocal bootstrapping variance of 
xp-* in the individual fields, Wij = l/a'^{x[^'^). This weight- 
ing scheme yields the minimum-variance estimate of the 
average Xi . 



Bootstrapping II. In a second bootstrapping stage, we ran- 
domly draw fields from the set of 52 GaBoDS fields (with 
putting back) and calculate a combined signal according 
to Eq. (|58l) for this bootstrap sample. We repeat this pro- 
cedure 2000 times. This estimates the PDF of statistical 
errors in the combined signal including cosmic variance 
and accounting for the adopted weighting scheme. In ad- 



dition to the combined signal of the aperture moments, 
we use this technique to also estimate the PDF of statis- 
tical errors of the bias parameters ^((^ap) and &(6'ap), Eqs. 
()50|l . which are computed from the weighted average of 
the aperture statistics. Based on the bootstrap samples 
we obtain the median, (asymmetric) 68% confidence in- 
tervals about the median and covariances of the errors for 
the final results. 



5. Results 

5.1. Aperture statistics 



The combined measurements for {N'^) 
can be found in Fig. El and Fig.[7| 



(TVAfap) and (Mfp) 



Galaxy clustering. Traditionally, galaxy clustering is stud- 
ied using the angular correlation function lu(0). For a 
comparison of our results for the two-point statistics of 
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Fig. 7. Top row panels and lower left panel: cross-correlation between aperture mass and aperture number count 
for the three different foreground samples FORE-I (solid boxes), FORE-II (open stars) and FORE-III (open crosses). 
The panels are subdivided; the lower panel shows the B-mode, upper panel is the E-mode of (NMnp). The curves are 
ACDM predictions assuming unbiased galaxies. Lower right: aperture mass dispersion, lower and upper panel are 
B-mode and E-mode, respectively. The solid line is a ACDM prediction. The solid lines in the B-mode panel are the 
E-mode prediction with positive and negative sign, which have been inserted for comparison. 



galaxy clustering with the literature it would be conve- 
nient to have uj{9) cleaned from the integral constraint. 
Already from early studies on galaxy clustering (e.g. 
Davis & Peebles 1983) up to recent studies (Norberg et 
al. 2001; Zehavi et al. 2002) it is known that oj{0) is close 
to a power law over a wide range, thus 

u;{0) = (^^) , (59) 

where A^^ and 6 are the clustering amplitude at 1' and the 
clustering power law index, respectively. Assuming this 
power law, we can calculate the aperture number count 
dispersion for our polynomial aperture filter: 

{N^)i9) = j\xxu;{ex)T+ix) = f{S)A^ (^^^ , (60) 

with f{d) being defined in (|A.3|I . Therefore, if uj{0) is a 
power law, the aperture number count dispersion is also 



a power law with the same slope but different amplitude. 
Since {N'^) is free of the integral constraint, we can recover 
uj{9) from (iV^) by fitting a power law and rescaling the 
best- fit amplitude by means of the function (|A.3p . For 
determining the statistical uncertainties of A^j, however, 
we have to take into account that it is not only a function 
of the amplitude of {N"^) but also a function of the slope 
S. 

We performed fits of l|Hnjl to the data taking into ac- 
count the covariances of (iV^) obtained by bootstrapping, 
see Fig. What is concluded for our foreground samples 
is summarised in Tabled 

The angular correlation of the galaxies in FORE-I - 
a sample roughly comparable to the foreground sample 
of Hoekstra et al. (2002) - has a slope slightly steeper 
than what is found in the sample of Hoekstra et al. (there 
S — 0.7 and A^^ = 0.115) and is smaller in amplitude for 
aperture radii larger than 0ap ~ 3'. This discrepancy in 
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Table 3. Amplitude and slope of the angular correlation, 
uj(9) = {9/l')~^ , in our foreground galaxy samples 
as inferred from (iV^)(0ap); X^/"- denotes the reduced 
{n =12 — 2) of the maximum-likelihood fit. 

galaxy sample 5 X^/'* 



FORE-I 0.086 ±0.006 0.89 ± 0.05 1.0 

FORE-II 0.044 ±0.003 0.99 ±0.04 0.6 

FORE-III 0.026 ±0.002 0.90 ± 0.04 1.8 
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Fig. 8. The measured two-point cosmic shear auto corre- 
lation in terms of ^± = (7t7t) ± (7x7x)- The dotted curve 
is a ACDM prediction based on ph{w) as in the sample 
BACK. 



j4(j and 5 is not as drastic as it may seem if one takes into 
account that the errors of A^^ and 6 are anti-correlated: a 
smaller A^ results in a steeper 5. Another issue that may 
play a role in this context is the fact that Hoekstra et al. 
use a different filter, R^, which is somewhat different from 
our i?-band filter. All in all we think that the measurement 
of u){9) for FORE-1 is consistent with the measurement of 
Hoekstra et al. 

Compared to the ACDM prediction of (iV^) for un- 
biased galaxies, which trace the dark matter distribu- 
tion, our measurements are clearly different, namely ex- 
ceeding the dark matter expectation on scales smaller 
than 0ap ~ 5', and falling slightly below the prediction for 
the largest aperture radii. This already suggests a scale- 
dependence of the bias factor. 

Dark matter clustering. The clustering of the total matter 
content as derived from the ellipticities of the background 
galaxies is expressed by the dispersion of the aperture 
mass, Fig. |3 We calculated this quantity for a range of 
different aperture radii from the cosmic shear two-point 
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Fig. 9. Plot of the measured mean tangential shear, 
(7t), about galaxies in our foreground samples FORE-1 
(filled boxes), FORE-11 (open stars) and FORE-111 (open 
crosses) as a function of angular separation. The lines are 
ACDM predictions assuming unbiased galaxies, Pi{w) and 
Ph{w) as estimated in our samples: for FORE-1 (solid), 
FORE-II (dashed) and FORE-III (dotted). 



correlators, ^j-, which are shown in Fig. |H1 (rebinned for 
that plot). 

In all figures, the prediction for the adopted fiducial 
cosmological model and the estimated redshift distribu- 
tions in our galaxy samples is plotted. We conclude that 
this prediction is in good agreement with our measure- 
ments. Therefore the fiducial cosmology taken for the bias 
parameter calibration seems to be reasonable. 

Judging from the B-modes, (M^), in Fig. [7| which 
serve as an indicator for systematics in the PSF correc- 
tion, the PSF correction is ok. Over the whole range of 
aperture radii considered the B-modes are consistent with 
zero, maybe with a minor exception at about 0ap ~ 3'. 
See Hetterscheidt et al. (2006) for a detailed discussion on 
this issue. 

Correlation between galaxy and matter distribution. The 

cross-correlation between the A'^-maps and the Map-maps 
is plotted in Fig. □ Apart from 6'ap « 3' in FORE-II the 
B-modes of the signal are all consistent with zero. The 
cross-correlation has been worked out on the basis of the 
mean tangential shear about galaxies in the foreground 
samples. Results for the galaxy-galaxy lensing signal are 
depicted in Fig. El 

The data points (E-mode) on intermediate scales are 
below the theoretical prediction for (A^Afap) based on an 
unbiased galaxy population. This again indicates that ei- 
ther the bias factor or the correlation parameter or both 
differ from unity, hinting towards a population of galaxies 
that does not perfectly trace the (dark) matter distribu- 
tion. 
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Fig. 10. The linear stochastic bias parameters of galaxies in the samples FORE-I, FORE-II and FORE-III (left to 
right column); the bias factor, &, is upper, the correlation parameter, r, is in the lower row. The parameters have 
been calibrated assuming — 0.3 and JIa — 0.7 (see Fig. 0}. The effective comoving scale is based on the aperture 
radius and the mean redshift of the weight functions /11/3, Fig. [51 The bias parameters for a particular aperture radius 
are averages over different physical scales and redshifts fSect. IT!5|) . The shaded area denotes the average bias factor 
or correlation factor over all aperture radii between 0ap = 2' . . . 19'; the maximum-likelihood of this average and its 
statistical uncertainty are shown in numbers inside the panels. 



5.2. Galaxy bias parameters 

The final result of our work is displayed in Fig. EH The 
bias parameters calculated from the aperture statistics, 
Eqs. (|50|) . have been calibrated, and the aperture radii 
have been converted into a typical physical scale, i?, based 
on the mean redshift of the range over which the param- 
eters are averaged. As this redshift range stretches over 
about 40% — 50% (Icr) of the mean redshift (see Fig. 
ISJ, there is a relative uncertainty attached to the phys- 
ical range, R, which is of the same order; for instance for 
R — 6 h^^Mpc we have as resolution for the effective scale 
an S/i^^Mpc (see Sect. ESI- 

Over the range of (comoving) physical scales investi- 
gated, below about R < 10/i^^Mpc, the bias factor stays 
more or less constant, rising towards smaller and possi- 



bly also larger scales with a valley on intermediate scales, 
where b becomes slightly inconsistent with 6 = 1 at a 
68% confidence level; this implies a scale-dependence of 
the bias factor. As absolute minimum we obtain b^dn — 
0.78 ± 0.10, 0.74 ± 0.10, 0.78 ± 0.10 at roughly 6'ap « 10'. 
The position of the minimum is not well defined, however, 
due its width. In order to get an average value for the bias 
factor, we make a maximum likelihood fit assuming a con- 
stant bias over the range 2' < (?ap < 19' while taking into 
account the covariance between the errors, as estimated 
from the bootstrap samples, shown in Fig. 1111 This fit 
yields: b = 0.81 ±0.11, 0.79 ± 0.10, 0.81 ±0.11 for FORE-I, 
FORE-II and FORE-III, respectively. Therefore, over the 
selected range of scales, galaxies are anti-biased, i.e. less 
clustered than the dark matter. 
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Fig. 11. Correlations of the statistical errors of the bias factor (left panel) and correlation factor (right panel) of 
galaxy sample FORE-II as inferred from the bootstrap samples; the correlation matrices of samples FORE-I and 
FORE-III are virtually identical. The colour of one pixel in an intensity map denotes the correlation between the 
errors belonging to the two bins defined by the x- and y-axis; the key of the intensity map is on the right side. The 
numbers attached to the axis denote the aperture radii in arcmin corresponding to the individual bins in Fig. 1101 



The correlation factor, r, has a larger relative un- 
certainty than the bias factor, b, since it is based on 
two lensing quantities, (iVMap) and {M^^}, which are 
generally noisier than (A^^). Broadly speaking, the cor- 
relation of the galaxies to the (dark) matter distribu- 
tion is relatively high. A scale-dependence of the cor- 
relation factor is hard to determine due to the large 
uncertainties and the high correlation of neighbour- 
ing bins; it may be present in the sample FORE-I. 
Averaging the correlation factor over 2' < 6'ap < 19' 
yields f = 0.61 ± 0.16, 0.64 ± 0.16, 0.58 ± 0.19 (FORE-I to 
FORE-III) which reflects both the high correlation and 
the unfortunately still large error bars. Obviously, a much 
larger survey area is required to obtain better constraints. 
We are going to discuss our results in the following section. 

6. Discussion and conclusions 

Observationally, the galaxy-dark matter bias can be 
probed by means of various methods (see introduction). 
Gravitational lensing provides a promising new method in 
this respect. It is special because it allows for the first time 
to map the total matter content (mainly dark matter) with 
a minimum of assumptions and independent of the galaxy 
distribution. Such a map can be compared to the distri- 
bution of galaxies, or particular types of galaxies, in order 
to investigate the galaxy bias. In particular, correlations 
between galaxy and dark matter density become directly 
visible. For working out the galaxy-dark matter bias, older 
methods rely on assumptions regarding the growth of dark 
matter density perturbations, the peculiar velocities of 
galaxies and their correlation to the dark matter density. 
Moreover, they often only allow one to measure the bias on 
large (linear) scales, > Sh^^Mpc, whereas the non-linear 



regime is also accessible with lensing. However, gravita- 
tional lensing has the disadvantage that it is not equally 
sensitive at all redshifts. The cosmic shear signal is most 
sensitive to matter fluctuations roughly half-way between 
z = and the mean redshift of the background. This de- 
fines a natural best-suited regime for the method at a 
redshift of about z « 0.5, often even slightly lower, con- 
sidering the depth of current galaxy surveys. It is expected 
that the most sensitive regime will be shifted towards 
higher redshifts by future space-based lensing surveys. 
Furthermore, lensing observables are quite noisy so that 
large survey areas are required for a good signal-to-noise. 
Impressively large surveys with instruments such as the 
CFHT (CFHT-Legacy-Survey, CFHTLS), the VST (Kilo- 
Square-Degree-Survey, KIDS), Pan-STARRS, or SNAP 
are either ongoing or about to start within the next years, 
providing us with plenty of high signal-to-noise informa- 
tion on dark matter and galaxy clustering. 

In this paper, we employed aperture statistics to quan- 
tify the relation between the dark matter and galaxy 
density. We tested the evaluation software against Monte 
Carlo simulated WFI fields, assuming an unbiased galaxy 
population, and found that the software is working to 
at least a few percent accuracy (Simon 2005). The data 
used is the GaBoDS with restriction to galaxies brighter 
than 24 mag in the i?-band; this allowed us to estimate 
the redshift distribution of the galaxies on the basis of 
three COMBO-17 fields (A901, AXAF/CDFS and Sll) 
for which photometric redshifts in < z < 1.4 are avail- 
able. For all the other fields, only i?-band magnitudes can 
be used to select galaxies. For this selection, we defined 
foreground galaxy samples by choosing galaxies from three 
i?-band magnitude bins that have increasingly fainter me- 
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dian magnitudes. The sample FORE-I is comparable to 
the foreground selection in Hoekstra et al. (2002) who ap- 
plied the same technique as we are using here. By means 
of the photometric redshifts of the COMBO- 17 fields we 
can translate a GaBoDS R-band magnitude interval into 
a redshift distribution. The fainter the bin, the broader 
the redshift distribution, while the mean redshift moves to 
larger values. Therefore, only FORE-I has a rather sharp 
peak in redshift, while FORE-III stretches between red- 
shifts of about z = 0.1 and z w 0.9. Hence, FORE-II and 
FORE-III are averages over a relatively wide range of red- 
shifts. In order to get narrower distributions in redshifts 
with the aim to reconstruct the redshift evolution of bi- 
asing, multi-colour lensing surveys are required. Cosmic 
variance is the main uncertainty in the estimated red- 
shift distribution. Based on the field-to-field variance of 
the photometric redshift distributions we estimate that 
this uncertainty translates into a Icr-uncertainty of the 
bias parameters of 7%, except for the bias factor, 6, in 
FORE-I which has - 16%. 

The B-mode of the aperture statistics (A^Map) and 
(M^p) are used as an indicator for systematics in the PSF- 
corrected shapes of the background galaxies (Fig. [7| ; they 
cannot be produced by gravitational lensing and should 
therefore be pure noise. We find that the B-modes are con- 
sistent with zero, maybe leaving a small question mark at 
6'ap ~ 3'. Note that, at least in principle, physical effects 
like intrinsic alignments of the source galaxies can be a 
source of B-modes in (Affp) (on small scales), so that van- 
ishing B-modes are not the ultimate indicators of PSF sys- 
tematics. For (A^Afap), however, the only possible source 
of B-modes is a violation of a statistical parity-invariance 
(Schneider 2003). Therefore, (TVMap) should always be B- 
mode free which is clearly the case in our data. 

The fit of a theoretical (Mfp) constructed from our 
fiducial cosmology and redshift distribution of source 
galaxies to the measured (M^p) is an important test for 
the calibration of the bias parameter; (Affp) is indepen- 
dent of the galaxy bias. Our data points are consistent 
with the fiducial cosmological model and, therefore, we 
accept the fiducial cosmological model and the estimated 
redshift distributions as sufficiently accurate for our pur- 
poses. For fiducial cosmological models different from ours 
(but still flat, rim + f^A = 1, with negligible baryon den- 
sity, r ~ ilm/i, and h — 0.7, and erg oc ilj^°-^^) the cali- 
brated bias parameters in Fig. 1101 may be scaled up or 
down using Fig.^ In a related paper (Hetterscheidt et al. 
2006), we discuss in much more detail issues concerning 
the creation of source galaxy catalogues and their data 
quality, and we determine constraints on cosmological pa- 
rameters based on the GaBoDS data. There it can be seen 
that this cosmic shear analysis supports our adopted fidu- 
cial cosmology. An uncertainty in the fiducial cosmology 
adds an additional uncertainty to the galaxy bias cali- 
bration and therefore the inferred bias parameters. For a 
realistic relative error of 10% in fim, we estimate this er- 
ror to ^ 7% for the bias factor, b, and to ~ 2% for the 



correlation factor, r. Errors given in the following do not 
include calibration uncertainties. 

The result of the galaxy bias measurement is plotted 
in Fig. ^1 Overall, the galaxy bias factor and the cor- 
relation are close to an unbiased population of galaxies, 
i.e. 6=1 and r = 1. A possible scale-dependence is indi- 
cated for the bias factor which rises to 6 > 1 on scales 
below 6'ap ~ 4', falls below 6 = 1 on scales of 6'ap ~ 5' and 
possibly rises again on larger scales. An aperture radius 
of 6'ap — 4' corresponds to an effective comoving scale of 
R = 1.4,2.1,2.8/i-iMpc (FORE-I to FORE-III) with a 
relative uncertainty (Icr) of about 40%. The origin of this 
uncertainty is due to the fact that we are actually observ- 
ing averages of galaxy bias over some redshift (cosmolog- 
ical time) and scale as illustrated by Fig. [SJ the median 
redshifts for the bias are z = 0.28, 0.44, 0.59, respectively. 

The median redshift for the correlation parameters 
are slightly different from those of the bias factor, here 
z = 0.30,0.38,0.43, again with relative widths (Icr) of 
about 40%. Thus, the correlation parameters reflect val- 
ues typical for a slightly different, more recent cosmologi- 
cal time. This mismatch arises if the peak redshift of the 
lensing efficency, W{w), is displaced with respect to the 
peak redshift of the foreground sample, as can be seen by 
Eq. H55(l . An alignment could be achieved by choosing an 
appropiate background sample for every foreground sam- 
ple which was not possible in our case, because we did not 
allow background galaxies fainter than 24 mag. 

Going back to the observed scale-dependence of the 
bias factor, galaxies become anti-biased on intermediate 
scales; they are less strongly clustered than the matter. In 
our data, the minimum value of the bias factor is deter- 
mined to be bmin ^ 0.76. This kind of scale-dependence 
has also been detected by Pen et al. (2003) (VIRMOS- 
DESCART survey) and Hoekstra et al. (2002) (VIRMOS- 
DESCART and RCS) which both rely on weak gravita- 
tional lensing to probe galaxy bias. While Pen et al. use 
/-band luminosities to select galaxies, which results in a 
larger value for the minimum bias factor but at a similar 
scale of about R w 3/i"^Mpc {k = 2tt/ R ^ 2hMpc~^), 
the data and sample selection of Hoekstra et al. is rel- 
atively similar to our sample FORE-I; their value of 
bmin = O-'^'l-o 04 agreement (Icr) with our measure- 
ment, but the quoted scale of i? « l/i~^Mpc is differ- 
ent. However, as emphasised before, the position of the 
minimum is not well defined in our data. Considering the 
statistical errors one has to admit that the position of the 
bias minimum is not well determined also in the Pen et al. 
analysis (their Fig. 19). Hence, there is no contradiction 
between our data and that of the other authors. 

An anti-bias on the scales considered here and a char- 
acteristic "dip" in the functional form of the bias fac- 
tor is in concordance with recent numerical simulations 
of dark matter structure formation (Springel et al. 2005; 
Weinberg et al. 2004; Guzik & Seljak 2001; Pearce et al. 
2001; Yoshikawaet al. 2001; Somerville et al. 2001; Jenkins 
et al. 1998). The scale-dependence is due to the fact that 
the galaxy clustering is a power-law over a wide range of 
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scales, reflected by (iV^) in Fig. El while the dark matter 
clustering has different shape in CDM simulations and in 
the observations suggested by, for instance, (Mfp) in Fig. 

m 

For the linear correlation parameter, we observe as 
Hoekstra et al. (2002) and Pen et al. (2003) a high corre- 
lation between galaxy and matter distribution. Averaging 
the measurement of Hoekstra et al. over the range 
2' < ^ap < 19' yields roughly r « 0.8 which is consistent 
with our average (Icr). Our observed correlations between 
fluctuations in the galaxy number and mass density ap- 
pear to be a bit lower, though (Hoekstra, private com- 
munication). This could hint to an hitherto undiscovered 
systematic effect in our data. However, it should be kept 
in mind that the statistical errors in r are highly corre- 
lated and quite large so that this slightly lower value of r 
may be just a statistical fluke. The clear scale-dependence 
of the correlation parameter observed by Hoekstra et al. 
is not visible in our data, because this feature probably 
gets lost within the statistical uncertainties. 

The figures for the correlation parameter - r is smaller 
than unity with 68% confidence - show that the galaxies 
are either stochastically or non-lincarly biased, or a mix- 
ture of both. To understand what is meant here, imagine 
that 6g ~ the galaxy density contrast - and (5m - the dark 
matter density contrast (both smoothed) - are quite gen- 
erally related by 

= f{Sm) + eiSnd , (61) 

where / is some function and e a random variable (noise) , 
both satisfying 

(/(^m)) = {<Sm)) = mn)f{Sn.)) = (62) 

owing to the definition of the density contrast and noise 
(e is statistically independent of /). In the case of / be- 
ing linear and e = we have a linear and deterministic 
relation between Sg and (5ni, whereas e 7^ introduces a 
stochasticity between the density fields. The latter case is 
called stochastic bias. If 6g and were Gaussian random 
variables, / would be a linear function. A non-linear func- 
tion / yields what is commonly called a non-linear bias. 
The degeneracy between non-linearity and stochasticity 
arises because a decorrelation - indicated in the linear 
stochastic bias scheme by r<l - can be generated by both 
a non-linear / and a stochastic component e: 

^ {SgS„,) ^ {SmfjSm)) .gg. 

Discriminating between these two cases requires the 
additional measurement of the non-linear stochastic bias 
parameter (Yoshikawa et al. 2001; Dekel & Lahav 1999) 
or equivalent quantities. So far, these parameters have 
only been measured for the relative bias between popu- 
lations of galaxies (Wild et al. 2005). Therefore, it is un- 
known from observations whether this decorrelation be- 
tween galaxies and dark matter is mainly due to a ran- 
dom scatter or a non-linear relation. To resolve this, using 



approaches similar to ours where statistical moments of 
the joint PDF of matter and galaxies are measured, one 
needs to invoke higher-order statistics. As the currently 
ongoing research is working on the three-point statistics 
of the aperture mass and the aperture number count (e.g. 
Schneider & Watts 2005; Schneider, Kilbinger & Lombardi 
2005; Jarvis et al. 2004; Schneider & Lombardi 2003) we 
can expect to be capable of such a task quite soon. 

Within the uncertainties of our measurement we do 
not sec a difference in the biasing parameters between the 
three foreground bins. As the three different foreground 
bins represent different median redshifts of the galaxies, 
we conclude that on the scales considered the redshift de- 
pendence of the (averaged) linear bias for 0.3 S z < 0.7 
has to be smaller than about Ab < 0.2 and Ar < 0.4 (Icr) 
as crudely estimated from the Icr-errors of the average 
bias and correlation parameter; any larger bias evolution 
should have been detectable despite the relatively large er- 
ror bars. These figures are no serious constraints to cosmo- 
logical models for the bias evolution because all different 
numerical and analytic models predict evolution rates well 
below these limits in the redshift range covered here (cf. 
Magliocchetti et al. 2000). In a recent paper, Marinoni 
et al. (2005) measured the non-linear biasing function 
(Dekel & Lahav 1999) in the VIMOS-VLT Deep Survey 
between 0.7 < z < 1.5 and found that the bias evolution 
is marginal below z 0.8 and becomes more pronounced 
beyond that redshift. Empirically, they found the redshift 
dependence of the bias factor on a scale of 8 h~^Mpc being 
described by b{z) ^1 + (0.03 ± 0.01) (1 + z)3-3±o.6 ^i^^^ch 
means a change of b of roughly 5% between 0.3 < z < 0.7. 
This figure is in qualitative agreement with our observa- 
tion. 

The bias parameters, no matter whether linear 
stochastic or non-linear stochastic bias, are just conve- 
niently defined quantities for a comparison of random 
fields. They bear no obvious relation to the physics of 
galaxies. In the end, these measurements will need to be 
interpreted in terms of physical quantities like the halo 
occupation distribution (Bcrlind, Weinberg et al. 2003; 
Berhnd & Weinberg 2002; Peacock & Smith 2000) in or- 
der to learn more about the evolution and formation of 
galaxies in the environment of their parent dark matter 
haloes. 
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Appendix A: Angular clustering and the aperture 
number count dispersion 

We will briefly demonstrate in this section that the inte- 
gral transformation H39|) applied to a offset power law 



c 



(A.l) 



where C is a constant (the integral constraint), results in 
a power law aperture number count dispersion {N'^). 

For the polynomial aperture filter u used here, Eq. 
H49|l. one obtains for the transformation kernel T+ an an- 
alytical expression that can be found in Schneider et al. 
(2002). Using this kernel and Eq. l|ITT|l for uj{e) into ^ 
yields 



ap 



(A.2) 



where the function f{S) has been obtained by 
Mathematical, 

m ^ 42 



3+S 
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30 900 



(A.3) 



i5r(^) 
2+^" 4+^^ V^r(3 + |) 



30 (2 + (5)V^ 



2r(^) 
r(i+|) 



(2 + <5)"0F 



1508r(^) 

%/^r(5 + |) 



ii6or(^) + 
V^r(4 + |) 



1056r(^) 288r(i^) 



V^r(7- 



Therefore, (iV^) is insensitive to the offset in uj{9) and 
is a power law with the same slope as Lij{9). In the 
regime 5g [0.2,1.6], the somewhat bulky function f{S) 
can, within a few percent accuracy, be approximated by 



f{6) « 0.0051 (5" '^^ + 0.2769 6^-^^ + 0.2838 S 



1.25 



(A.4) 



which covers the commonly observed range of values for 
the power law index. 
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